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EJECTION OF INTERNAL MASS 
FROM SPINNING SPACECRAFT 
WITH COLLISION INTERACTIONS 
Franz Zach 

ABSTRACT* 

This report treats ejection of internally stored mass (e.g. an apogee motor 
casing) from a spinning spacecraft wherein a collision occurs. The problem is 
attacked by formulation of twoprocesses: first, an ejection process which includes 
the release mechanism, the separation trajectory of stored mass relative to 
spacecraft, and the determination of the collision point; and second a collision 
process wherein the impulse exchanged is computed based upon consideration of 
surface friction, mechanical deformations , and coefficient of restitution. A spe- 
cific case is treated where both direction and magnitude of spacecraft body rates 
following collision are determined as functions of the collision parameters. A 
knowledge of these body rates is required to perform directional despin of a 
spacecraft about its final axis of maximum moment of inertia. 


*The work for this paper was performed while the author was granted a National Research Council 
Postdoctoral Resident Research Associateship supported by the NASA/Goddard Space Flight 
Center. The author is very much indebted to Mr. W. Isley, Mr. E. Stengard, Mr. R. Bartlett and 
Mr. D. Endres for their contributions. 
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r 


distance from CMS to CMC 
brackets of the solar aspect sensor 
distance from CMM to CMC 

center of mass of spacecraft and internal mass together, before 

ejection of internal mass 

center of mass of internal mass alone 

center of mass of spacecraft alone 

angular momentum 

energy transferred from motion into deformation 
coefficient of restitution 

unity vectors in directions x, y or z respectively 

range of spring operation 

impulse 

moment of inertia matrix 
maximum moment of inertia 
motor 

mass of the internal mass alone 

mass of the spacecraft alone 

unity vector normal to collision surface 

point of collision on the spacecraft measured in the Xg , y s , z s set 
point of collision on the motor measures in the x M , y M , z u set 
distance vector 

distance vector from CMS to P 
distance vector from CMM to P 
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spacecraft 
collision vector 
solar aspect sensor 

spin mode before deployment of internal mass 
spin mode after deployment of internal mass 
time 
velocity 

velocity of CMM 

velocity of colliding point of motor 
velocity of CMS 

velocity of colliding point of spacecraft 

right handed coordinate system with center in CMC and axes along 
axes of principal moments of inertia 

right handed coordinate system with center in CMS and axes parallel 
to the x, y, z set 

right handed coordinate system with center in CMM and axes 
parallel to the x, y, z set 
coefficient of friction 
inertial coordinate system 

coordinate system with center in CMM and axes parallel to £, t }, £ 
angle of friction, = arc tan 
angle of S against -y 
vector of angular velocity 
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EJECTION OF INTERNAL MASS FROM SPINNING SPACECRAFT 


WITH COLLISION INTERACTIONS 


1. INTRODUCTION 

A spacecraft despin maneuver is often required following orbit injection 
and deployment of the apogee motor casing. This report treats ejection of such 
stored mass under conditions where the spin axis after apogee motor burnout is 
not the spacecraft axis of maximum inertia. This condition could result from 
a damping mechanism which is not compensated for by active nutation control. 

Under such circumstances, the spin axis will now be aligned along the new 
axis of maximum inertia which will not be the axis of symmetry, normally 
assumed for motor ejection. This can produce a collision process between the 
ejected mass and spacecraft, wherein the resulting body rates can dictate the 
final spin direction. The determination of body rates can be approached through 
definition of an Ejection Process and a Collision Process. In the Ejection 
Process, the release mechanism, departure trajectory and collision point are 
established. This is then followed by formulation of the Collision Process wherein 
impulse exchange is computed in terms of surface friction effects, mechanical 
deformations, and coefficient of restitution. Based upon the impulse exchange, 
direction and magnitude of spacecraft body rates are determined. Some despin 
mechanisms, such as Yo-Yo devices , require a specified initial spin direction 
in order to remove angular momentum. If the axis of deployment symmetry 
becomes the axis of maximum inertia following motion ejection, then final spin 
direction will be determined by direction and magnitude of body rates after 
collision and damping characteristics of the spacecraft structure. 
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It is therefore apparent that estimation of body rates after collision can 
dictate ultimate despin capability of the spacecraft. Part 2 develops equations 
for the Ejection Process. The release mechanism is modeled in Section 2.1 as 
an impulsive device, which essentially idealizes the acceleration as a pure 
addition of velocity. 

In Section 2.2 the release mechanism is defined to represent time dependent 
acceleration. In Section 2.3, Collision point and closing velocity are determined. 
Part 3 treats the collision process. Section 3.1 considers a representative static 

—4 

structural load test. The collision impulse vector (S) is computed as a result 

of much test data. In Section 3.2 a mathematical model is developed for calcu- 

— ♦ 

lation of S. This Section also covers determination of body rates due to the 
collision process which are evaluated as functions of coefficient of restitution 
and surface friction. 

2. The Ejection Process — Trajectory and Collision Point 

2.1 Release Mechanism - Impulsive 

Figure 1 shows an example of the geometric configuration for the space- 
craft and stored mass which is represented herein as an apogee motor casing. 

It shall be assumed, that the x-axis is the axis of maximum moment of inertia 
and that spin is stabilized about this axis. This fact is indicated by in Fig. 1. 
It shall be assumed further that the apogee motor is released during this rotation 
and accelerated by springs in the z-direction. The springs usually are effective 
only over a small distance and are located rotational- symmetric about the 
z-axis. It is assumed herein that the release mechanism involves no tip-off 
errors. This means that there are no reaction torques between spacecraft and 
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apogee motor due to separation, and implies that no change of angular velocity 
is produced during separation. Therefore, both spacecraft following separation 
and apogee motor casing will have the same Z as the original composite space- 
craft. The fact that the springs usually work only over a small distance results 
in a simplification of the problem: The motion of spacecraft and apogee motor 
can be calculated under the assumption that the springs Impart a relative 
velocity in z-dlrectlon to spacecraft and apogee motor in an impulsive manner. 

Figure 2 shows the situatior. immediately after release. The relative 
velocity in z -direction is given by the experimental data for the spring force and 
shall be called v ff : 


v 

r i 



'S* 


(i) 


Because there are no other forces active, the law of conservation of impulse 
yields 


"Hi I v mJ " ms I v s«l = 0. (2) 

which together with (1) allows calculation of | v Mi | and I v Sf | . v My and v Sy 
are obtained as follows: upon separation, the CMM rotates through displacement 
b about CMC with the angular velocity oj % , CMS rotates about CMC with co t 
through displacement a (see Figs. 1 and 2). Because there are no impulses in 
y-direction due to separation, the motor will move with velocity: 

v „ = - a; x b (3) 
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and the spacecraft will move with velocity 


As some time t after deployment, relative position and motion in an inertial 
coordinate system will be as shown in Fig. 3 . The displacements are given by: 

’i (o = 

Ci(t) = C ll (t).C i i.)J 

and 

y (t) = - i* (t) sin 0)t + tj* ( t) cosa > 
z (t) = i* (t) cos 'ot + 77* ( t ) sin cut 



y(t) and z(t) give the motion of the apogee motor with respect to the CMS. 

Based on Figs. ( 1 ) through ( 3 ), displacements in the 77, £ coordinate 
systems become: 

T?„(t) = v My t = I oj x | bt, 

C M (t) = v Ml t+b, I 

V s (t) = v Sy t = - |oj x | at , 

m 

= V S. 1 - 3 = “ V M,-- 1 " a - 

m s J 


7?*(t) = (a + b) |«J t "I 

and / m \ > (8) 

4 ,I(t) = v m* + -~j t + a + b = v rt t + a + b.J 

Entering Eq. (8) into Eq. (6) now gives the time history of the CMM in the y-z 
coordinate system. Because the same u j g is valid for spacecraft and motor 
after deployment, Eq. (6) produces the trajectory of every part of the motor 
with respect to the spacecraft, the initial condition being given by Fig. 2 . The 
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axes of rotational symmetry (aligned along the z-axis in Fig. 2) of spacecraft 
and apogee motor always remain parallel. 

Figure 4 shows trajectories based upon the following values: 


a = 3.2", b = 26.24" ^ 

a> x = - 10 rpm to - 150 rpm> 


0 ) 


v rz =41.9 in/sec . J 

Calculations with varying co x and v f2 = o produce trajectories which fall 
on top of the case where co x = -150 rmp in Fig. 4. This shows the spring to 
have little effect on the trajectory foraJ x greater than 92 rpm. 


2.2 Release Mechanism - Time Dependent Acceleration 

The previous section considered the release mechanism to be impulsive. 

For time dependent acceleration consider Fig. 5 and assume at first deployment 
without a spring. The deployment is assumed to be in the Instant when the z-axis 
passes through. the £' -axis. The 17 ' - £' system has its origin in CMS, the 
77' axis is not rotating and points in the z-direction in the instant of release. 
Defining r(t) in terms of Eq. ( 8 ) gives 


r(t) 



(10) 


and 


r(t) = (a + b) (1 + t 2 )- 3/2 


( 11 ) 


For the additional acceleration caused by the spring the force F is: 

F = K [r - (a + b)] 

where K is the spring constant. 
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The acceleration of the apogee motor due to the spring is F/m M and the 
acceleration of the spacecraft is F/m s . 

The sum of these accelerations gives if (t) due to the spring. The validity 
of this calculation is over the range f of spring- operation. 

Entering the additional acceleration into Eq. (11) gives 


r(t) = (a + b) w* (1 + a)*t)~ 3/2 + K — - - [r - (a+b)] 

m M m s 

for a + b ^ r - a + b + f . 


Consider a specific case having the following parameter values: 


i. { = -92 rpm ^ 

K = 2.81 x 10 6 slg s~ 2 
f = 1" 

a + b = 29.44" J 


(12) 


(13) 


Calculations of the two extremes, i.e. 

(a) no spring at all 

(b) spring action executed impulsively 

show approximately no difference in the trajectories. Eq. (12) gains importance 
for higher spring forces and lower |ojJ . 


2.3 Collision Point and Closing Velocity 

The trajectory of the CMM has been developed in Sections 2.1 and 2.2. The 
determination of initial collision point must be based upon the surface configu- 
ration of motor and spacecraft. 

In Fig. 6 0 i (x, y, z) represents any selected point on the motor surface and 
P i (x, y, z) a corresponding point of intersection with the spacecraft surface as 
determined by transformation of the CMM trajectory into Q i (x, y, z) The 
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initial collision is found by solving the minimum time problem for all per- 
missible Q . (x, y, z ) mapped into the set of all points P i (x, y, z). 

A geometrical solution is most often the simplest means of finding the 
initial collision point where surface envelopes are matched in a graphical 
maimer. 

The solution for the initial collision point shall be designated P = P(x gl y z g ) 
on the spacecraft and Q (x Mi y Mi z^) on the motor. 

Xg, y s , z s coordinates are measured with respect to CMS. The x M , y M , z M are 
coordinates measured with respect to CMM* Consider a configuration which 
yields a collision point characterized by P = (-3.9,- 22.5, 32.9) and Q = (-4.8, -16.1, 
-6.3). The time to collision at point P is obtained by solving Eqs. (6) and (8) 
for t with 


y(t) = Ay = y g - y„ 


and 


z(t) = Az = z s - (a + b + z M ). 


The closing velocity components can be obtained by differentiation of Eqs. 
(6) and (8) with respect to time. Using the following values: 

y s =- 22 - 5 " 

z P = 32.9" 

z m = - 6-3" 

the closing velocity components are found to be: 
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v =-171 in/sec 
yr 

and 

v zr = 164 in/sec. 

3. Collision Process 

Once it has been established that a collision takes place using methods 
described in Section 2.3, it becomes necessary to formulate a dynamics 
model for the collision process. 

It is important to calculate the impulse exchange between internal mass 
and spacecraft whicn gives the initial condition of translatory and rotational 
velocity after collision. 

—4 

For this purpose a good approach is to assume a collision vector S 
(Ref. 1) which is the time-integral over the forces in action during collision: 

S = f Fdt. (14) 

F as a function of time is given by the geometry at the collision point, by 
the amount of the relative velocity, by the masses of the colliding bodies and 
the structural properties of the two bodies in the area which is effected by the 
collision. 

For the purpose of evaluating a specific collision process, it is assumed 
that the motor and spacecraft configurations at the instant of collision are as 
represented in Fig. 7. The impact geometry shows initial contact with a 
bracket (1^ ) which supports a solar aspect sensor at 12° angle measured from 
minus y. The second impsct takes place at bracket (B 2 ) at 30° angle. 
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3.1 Collision Vector Determined From Static Structural Data 


Assume that static measurements have been performed for simultaneous 

loading of two brackets of the type shown in Fig. 7 and that deformation 

behavior is as represented in Fig. 8. From Fig. 7 geometry, bracket B 2 will 

not be contacted until bracket Bi has deformed .625". Let t = 0 denote the 

time when Bj deformation is initiated, t = tj the time when B 2 deformation is 

initiated, and t = t the end of deformation. The time t is also the moment 
2 2 

when the entire kinetic energy, present at t = 0, is absorbed by the deformed 
material. Fig. 9 gives the deformation diagram for both brackets with the time 
history outlined above, where t 2 occurs when the area under the deflection 
curve for B t + B 2 equals to the initial kinetic energy. 

Using the value of Vyr taken from Section 2.3, the kinetic energy to be 
absorbed by the structure is: 


T = 4618 lb in. ( 15 ) 

The kinetic energy absorbed at t can be taken from Fig. 9 as 230 lb-in. 
Therefore, 4388 lb- in., are to be absorbed when both brackets are deformed 
simultaneously. It should be noted that the deformation of the first bracket is 
already in a more advanced phase than the deformation of the second bracket. 

Considering Eq. (14) for S the direction and magnitude of F has to be 
given at all times. The magnitude is given by Fig. 9, the direction by the 
surface at the collision point, and the friction coefficient This latter term 
gives the maximum angle of deviation of actual force acting between the two bodies 
from normal to surface at collision interaction. There is always some 
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uncertainty in predicting the surface configuration due to (1) reproducibility of 
test data and (2) proper modeling of the structural dynamics. However, in' 
this case a good assumption can be made by knowing that the surface of the 
motor at the collision point is highly rigid whereas the brackets are deformed 
very easily. The collision surface therefore is given primarily by the shape 
of the motor. Considering Fig. 7 shows a deviation of 12° from a plane normal 
to y during deformation of For deformation of brackets Bj and B 2 both 
are considered separately, therefore giving the same 12° as above for B x and 
30° for B 2 . Fj will be called the force acting on Bj, F 2 the force acting on 

B 2- 

According to Eq. (14) we receive with 

ds 


yr dt 


= v(s) 


where s is the deflection taken from Fig. 9: 


S = J F dt = J* F —L 
v yr 


(16) 


where v r can be calculated as follows: the kinetic energy at t = 0 is given by 


T (0) = - m„ v 2 (0) 

y v ' 9 M yr v ' 


(17) 


and will be reduced due to the area E under the curve in Fig. 9: 


T(0) -T(t) = E.. 


(18) 
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With 


T (t) 
y 



( 19 ) 


we receive 


V yr(0 = 



(20) 


which is given in Fig. 10. 

Determination of the direction of F in every moment involves the following 
considerations: with very high friction the force would be aligned along y, 
because this is the direction of impact. With small friction deviation of the 
direction of the force from normal to collision surface cannot be different 
more than given by/z. 

The normal to the impact surface at Bj is 12° and at B 2 is 30° off from 
y-direction. If we assume that the bracket is made of aluminum, the motor 
of steel, and the contact surface is lubricated, ^ will range from .1 to .2. 

Then, F 2 will deviate from the y-direction by 24.3° for ^ = .1, and 18.7° 
for m = .2. 

Eq. (16) can be written in the form 


S = 


f ,(t i ) F f s(t 2 ) p f s(t 2^F, 

J 7 ds+ J 7 ds= J T ds + 

*(°) s(tj) S (0) 


r ,(t 2> f 

i — ds (21) 

J v 


s(tj) 


since F 2 is i 0 only for t t < t < t ^ Figure 11 now shows F/v for B j and 
for B 2 , both as a function of the deflection, s. 
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The area under the graphs in Fig. 11 now gives the absolute value of S 

for I?! and B 2 . However, it must be noted that this is only valid for an inelastic 

collision process, where no forces are acting between the two bodies after the 

maximum deformation of the structure is obtained. 

(The direction of S for B 1 and B 2 is given by the normal to the collision 

surface at Bj or B 2 respectively and by/u. ). 

The area under the graphs in Fig. 11 can be determined graphically from 

s = 0 to s = 2.1". Fj/v and F 2 /v reach infinity at s = 2.17"; therefore, from 

s = 2.1" to s = 2.17" an analytic form of integration must be applied. 

F and F constitute the numerator for F /v and F /v and only the interval 
12 12 

2.1" - s - 2.17" is considered F and F can be treated as constants for the 

1 2 

purpose of integration. 

For v(s), Eq. (20) can be applied where the term -2E d /m M can be expressed 
in terms of s. This can be done by graphic approximation of E d in Fig. 10: 
from s = 1.75" to 2.5" E can be approximated by the linear expression: 

E d = — ( k 2 s " k i>- (22) 

m M 

Evaluating k, and k„ from E J in Fig. 10 and substituting in Eq. (20) yields: 

12 d 

v — v 0 /2.36 - 1.09 s. (23) 

Integration in Fig. 11 now results in: 

|S, | = 37.6 lb-sec "I 

and - r (24) 

I S 2 | = 16.7 lb-sec. J 
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We are interested mainly in the further motion of the spacecraft, namely 
^immediately after collision. For this purpose the change of angular momen 
turn shall be calculated based upon 


AD = 



r x F dt 


(25) 


where 0 - t - 1 2 is the time interval during which the collision happens, and 
generally is assumed tobe very short, r is a vector leading from the CMS to a 

—4 

point where F is acting. Here, two points of collision must be considered with 

4 * —4 

r t leading to Bj and r 2 leading to B 2 , where ^ and r 2 are assumed to remain 
constant during 0 - t - t 2 . Therefore, Eq. (25) can be written as: 


AD 


t j t 

= r x x J Fjdt + fjxj* 


Fj dt = Tj x Sj + r 2 x S 2 


Using Fig. 7, f x , r 2 , Sj and S 2 can be determined as shown in Fig. 12. The 
collision point can be represented by: 


r, = - 3.5 e x - 21.9 e y + 32.9 e, 


and 


r 2 =-8.4e x -20 e y 


+ 32.9 e z J 


(26) 


Using Fig. 5 it can be determined that a collision takes place when the tra- 
jectory of the motor forms an angle of 45° with the z-axis. Therefore, the 
following expression for components of the collision vector is valid: 


S 


1 Z 



for i = 1, 2 . 


(27) 
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T C/i T C/3 T C/3 3 IW 


Now the following expressions can be written with p = arctan p 
, = -37.6 sin (12° - p) e x - 37.6 cos (12° - />) e y + 37.6 for 0^12° 
j = - 37.6 e y + 37.6 e f for p Z 12°, 

2 = - 16.7 sin (30 ° - p) e x - 16.7 cos (30° - p) e y + 16.7 e z for 0^p^30°M 28 ) 


id 


- 16.7 e y + 


16.7 e for p - 30° , 


j 


We are mainly interested in changes of angular velocity A D § of the space- 
craft about its z-axis. Combination of Eqs. (25), (26), and (28) yields: 



l . i 


This leads to 


AD i# = 1161.5 sin p - 164 cos p for 0 £ p - 12°, 


y <**) (29) 


AD fi = 90 + 339 sin p - 80 cos p for 12° - p - 30° 


(c) 


and 


AD ig = 190 for p ± 30°. 


J 


(d) 


The above mentioned uncertainties in surface configuration can be put into 
Eqs. (25) through (29). For example, the uncertainty in p can be considered 
by changing p in Eq. (29) to smaller and greater values starting from the 
most likely value of p . It can be calculated from Eq. (29b) that 
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and 


AD I# - 0 for n = t gp - . 42 
AD fi 2 0 for ii = tgp * . 142 J 



(30) 


If a positive value of u after collision is desired, n = .142 gives the 
lower limit for friction between the colliding surfaces. 

As pointed out earlier, the above calculations of | S | by integration of the 
areas in Fig. 11 are made under the assumption of inelastic collision. To 
consider different degrees of elasticity a different approach is applied as 
follows. 


3.2 Collision Process Based Upon Coefficient of Restitution 

In this section we apply another approach for calculating the collision 
vectors Sj and S 3 . For simplification, Sj and S 2 are combined to one vector 
S. S can be calculated directly based upon the dynamics of the collision when 
a coefficient of restitution is assumed. The coefficient of restitution gives the 
degree of elasticity involved in a collision process. The coefficient of restitu- 
tion e (Ref. 2) is defined as 


e = 


V S ~ V M 
V M- V S 


(31) 


where v M and v s are the velocities normal to the collision surface before 
collision of motor or spacecraft, respectively, and v^ and v^ are these values 
after collision. Generally, all prime terms represent values after collision. 

All velocities have to be taken at the point of collision and are different in 
general from the velocities of CMS and CMM. 
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e = 0 denotes a completely Inelastic collision, 
e = 1 a completely elastic collision. 

In order to receive the change in motion due to collision, we relate ';he inputs 
and the angular momentum of both the spacecraft and the apogee motor to the 
collision vector. 


V - Is = S 

V - I M = -3 

V -^s = r s * S 

= r M X 


•\ 


> 




(32) 


r ( is the distance vector from CMS to the point of collision P, r M is the di^lance 
vector from CMM to P. Consider only the change of motion due to collision and 
set f s = 0 , D m = D g =0. 

Eq. (32) can be written as: 


and 


m s v s r ^ ' 

— • 

"Vi &u - v m> =" S ’ 

[l s^ w s = r s x S 


(33) 


The velocity immediately before collision is v Mp = v,, , and v s = v SP «= 0. 
The velocities immediately after collision are: 
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and 


(34) 


v s ' P = Vs + [Is]" 1 [ r s x S] X r s 



= v + 

M 




-1 


[7 b x (-s)i x ; m . 


With n denoting the unity vector normal to the collision surface, Eq. (31) can be 
written as follows: 




iti s ] 


-l 


(r xS)]x r' 


- v 


- [[I ]' 

M 


( x(S))] x 


r } • {n} = 


(v } 

M 


{n} 


(35) 


(x denotes a cross product, * a dot product). 

Instead of assuming two different collision points (IJ and P 2 in Fig. 12) we 
assume the collision impulse acting at a point P between P t and P 2 . S is now 
determined by its magnitude | S | and angle a against the minus y-axis. The 
amount | S| , can be calculated from Eq. (35) by assuming e. It is found that for 
assuming full inelasticity | S| based on this calculation is equal to I Sj + S 2 1 as 
calculated in Section 3.1. The angle a is determined herein on the basis of an 
inelastic collision process, where it has been tacitly assumed that S produces the 
same resulting D si as in Section 3.1. This condition results in Fig. 13 v/here 
Eq. (29) is applied and compared with 

A D = T s x S, (36) 


where 
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(37) 


r s = - 3.9 e x - 22. 5 e y + 32.9 e x 


is assumed as distance from CMS to P. 

Based upon & taken from Fig. 13 and the condition according to Eq. (27) 
S can be expressed by: 


S O • 

- - sin <r e - 
72 


S S 

ri cos a *7! *■ 


where 


S = | S • 


> (38) 




The structure Is deformed continuously during collision. It is reasonable to 
assume that deformation Is such that n , (which is necessary for application of 
Eq. (31)) lies approximately in the direction of S. 

Therefore, Eq. (38) oa.i be written as 


S = Sn 


(30) 


and 


n 


sin <7 

V2 * 


V2 Cy 



(40) 


This flexibility in considering structural deformation at the point of collision is 
an advantage of this approach compared to that used in section 3.1. 

Eqs. (33) and (39) now are inserted into Eq. (35) in order to receive the 
following equation for ! S | : 
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s = 


(e + 1) n * Vj 


|j^n + [ [Ig ]" 1 (r g x n)] x r g - -^-n - [ [l„] -1 [r M x (- n)] ]x r„j 


(41) 


{n> 


S now can be inserted into Eq. (33) and the motion of spacecraft and motor after 
collision can be calculated. 

For this purpose the following special values for masses and moments of 
inertia are chosen: 

a. for the spacecraft and apogee-motor combined 


M = 1125 lbs. 

C 

I = 115.6 slg ft 2 

I =111.9 slg ft 2 

I C8I =105.9 slg ft 2 

b. for the spacecraft alone 


M = 1003 lbs 

s 

I = 87.7 slg ft 2 

Sxx 

I Syy =84.0 slg ft 2 
I. = 101.17 slg ft 2 

c. for the apogee-motor casing alone 
M m = 122 lbs. 

I = 7.35 slg ft 2 

Mxx 

I u =7.35 slg ft 2 

I = 4.33 Slg ft 2 

Mzx 

r g is given by Eq. (37) 
r = -4.8 e -16.1 e + 6.3 e 

M x y z 
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Figs. 14a and 14b show the change in u> x , w y and <w f for the spacecraft as 
^unction of cr and e. The value of cr t where changes its sign is called <^ criticil 
The corresponding value for n can be found to be M critic ,i = .141. Lower 
frictions than a , give negative w contributions to the motions of the 

'critical ° ° * 

spacecraft after deployment of the apogee motor. 

4. Discussion of Results 

Equations developed in Sections 2 and 3 can be used to establish ejection 
and collision processes if the following parameters can be specified: 

1. Moments of inertia for spacecraft and ejected body 

2. Mass of spacecraft and ejected body 

3. Geometrical configuration for C.M. locations in spacecraft and ejected 
body 

4. Surface friction properties at impact point 

5. Structural elastic properties at impact point (expressed in terms of 
coefficient of restitution) 

6. Release mechanism behavior 

7. Velocities and angular rates prior to release 

8. Impact geometry and deformation time- history 

Usually, parameters (1), (2), (3) and (6) are readily available for a given 
spacecraft. Parameter (7) can be established by in-flight measurement. The 
uncertainties associated with these items are considered herein to be small 
and would have no significant influence upon the impact study. However, param- 
eters (4), (5) and (8) are generally not known to such precision and must be 
treated within a bound of uncertainty in order to predict limits of impact be- 
havior. For example, both surface friction and elastic properties are normally 
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based upon practical tests wherein reproducibility of data is of prime concern. 

The uncertainties associated with these parameters can be used to determine 
a range of collision vectors, S, which then permits a prediction of the upper 
and lower bounds on resulting body rates following impact. 

Consider now the special case described in Sections 2 and 3 where param- 
eter (8) is based upon static test data. The influence of parameters (4) and (5) is 
shown in Figures 14a and 14b. The friction coefficient for a lubricated surface 
of the selected materials is estimated to lie at a value between 0.1 and 0.15. 

The angle sigma can be used to set upper and lower bounds of uncertainty on 
the friction coefficient. The upper bound on sigma is 17.5 degrees, which 
results from the consideration that friction coefficient becomes zero at this 
angle. The lower bound on sigma is 0 degrees and represents the point where 
friction coefficient exceeds a value of 0.577. This limiting value is based upon 
the angles shown in Figure 7, wherein tan 30° equals 0.577. A range of values for 
coefficient of restitution, e, from 0 (inelastic case) to 0.5 was considered for the 
structural elastic behavior. Higher values for e are not likely to occur for the 
subject case. Figure 14b shows that e has no effect upon determination of 
direction for oj . For the most likely condition of inelastic collision, u t is 
found to range from +0.12 radians per second (6.87 degrees per second) at friction 
factors of 0.577 and above to - 0.10 radians per second (-5.73 degrees per second) 
at zero friction. Where no lubrication is applied to the contact surfaces, one can 
expect the friction coefficient for bare metal surfaces to be unity or higher. If 
the surfaces are painted, one could expect even higher friction factors. Either of 
these latter two conditions would produce an u of about +0.1 radians per second. 
It is of interest to note that the lubricated condition will produce an &> t having 
bounds of +0.01 radians per second and -0.028 radians per second. 


21 


5. Conclusions 


A semi-analytic approach has been developed herein for determining the 
residual spacecraft body rates produced by impact with an ejected mass. Two 
processes are involved. The first phase describes the mass ejection, which is 
handled analytically. The second phase is the collision process which requires 
formulation of an empirical deformation model which is combined with an 
analytic model for the collision vector. 

A specific case, examined using this mathematical approach, revealed that 
the residual spacecraft z axis rate can be bounded by coefficient of restitution 
and surface friction properties. The coefficient of restitution influences only 
the magnitude of the residual rate whereas the direction is given by surface 
friction. 
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Fig. 14b. Changes of cj t as a function of cr (and therefore as a function of m 

according to Fig. 12). The parameter is the coefficient of restitution 
which gives the degree of elasticity involved. 





z/in 




28 


i >, -him Ill' i llllillllllilllllHIHHWW*-* tinwwwmwwiluil'l" Ill* 





29 


Figure 6. Collision point mapping. 
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Figure 8. Static load test for deformation of both brackets Bj and Bj, both starting 
being deformed at the same time. 
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Figure 11. Ratio of F ( (i = 1, 2) and closing velocity versus deflection s for 
calculation of the collision vector S. 
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Figure 13. Angle a of equivalent collision vector S against negative y-axis compared with 

/O for ^ and 5 2 . 
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Figure 14a. Changes of cu x and as a function of cr (and there- 
fore as a function of /i. according to Fig. 12). The parameter is 
the coefficient of restitution which gives the degree of elasticity 
involved. 
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Figure 14b. Changes of co* as a function of cr (and therefore as a function of fJ. according to 
Fig. 12). The parameter is the coefficient of restitution which gives the degree of elasticity 
involved. 
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